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In a previous article concerning image distortion in non-perturbative gravitational lensing theory 
we described how to introduce shape and distortion parameters for small sources. We also showed 
how they could be expressed in terms of the scalar products of the geodesic deviation vectors of 
the source's pencil of rays in the past lightcone of an observer. In the present work we give an 
alternative approach to the description of the shape and distortion parameters and their evolution 
along the null geodesic from the source to the observer, but now in terms of the optical scalars (the 
convergence and shear of null vector field of the observer's lightcone) and the associated optical 
equations, which relate the optical scalars to the curvature of the spacetime. 



I. INTRODUCTION 



The distortion of images caused by the gravitational field of a massive deflector (or "lens" ) is very well understood 
in the case of weak fields and thin lenses, where the gravitational field can be treated via the linearsuperposition 
of the fields of point-like masses. In this case, imagesappear multiplied, magnified, or distorted, depending on the 
alignmentof the source, the deflector and the observer. We assume the reader isfamiliar with the standard, thin-lens 
theory of gravitational lensing asis found in [jjj. 

The standard theory of gravitational lensing provides the necessary tools for the study of the distortion of images 
via the linear mapping Y — AX where A is the Jacobian of the lens map, 

y = F(x*)=x + ^a(x) (1) 

i.e., where A = dF/dx and where x is the angular position of the image (as seen on the observer's celestial sphere) 
and y is the angular position of a point source in the background (or unlensed) spacetime. Di and D s are the distances 
of the observer from the lens and source and Di s is the distance between the lens and source. We denote by d(x) 
the deflection angle of the lightray in the lens plane. Note that both x and y can be considered, respectively, as the 
rescaled position coordinates, £ = xD\ of the image and ff = yD s of the source, in the lens and source planes, as 
represented in Fig. El When no lens is present, the lens equation reduces to y — it and A = I, the identity matrix. If 
there is a small, but extended source, we denote a central point of the source by yo and any point on its boundary by 
y. The central point is imaged at xq, while the image of the boundary point is at x. The displacement of the image 
boundary from its center, Ax = x — xq = X and that of the source, Ay = y — tjq = Y, are related by the Jacobian of 
the lens map, i.e., by 



The map is frequently expressed as 



Y = AX. (2) 



(3) 




where the quantities k and 7 = VjF+tI are interpreted, respectively, as the convergence and shear of the image 
wih respect to the unlensed source. [Later we will give a more precise meaning to these quantities in terms of the 
convergence and shear of the null geodesic congruence defined by the observer's past lightcone.] The inverse of A, 
i.e., A" , is referred to as the magnification matrix and "carries" the source's shape into the image's shape, via 

X = A~ l Y (4) 
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We are here interested in obtaining an analog of this approach to distortion in a generic case that applies regardless 
of the strength of the gravitational field and without reference to thin lenses. We refer to this as a non-perturbative 
approach to image distortion. 

In this paper, we consider small elliptical sources, which can be completely described by three shape parameters: 
area, semiaxes ratio and semiaxes orientation relative to a fixed direction. They are dealt with by means of connecting 
vectors in the null geodesic congruence that forms the past lightcone of the observer. The source and image descriptions 
in terms of Jacobi fields are developed in full detail in our preceding paper where we obtain expressions for the 
three parameters that measure the "total" distortion of the image with respect to the source: the image's solid angle, 
semiaxes ratio and orientation as compared to the source's parameters. 

In addition, we consider the distortion of the pencil of rays of the elliptical source as it travels towards the observer, 
where it eventually defines the image. This is an "infinitesimal" distortion; The distortion of the image is the result 
of the cumulative effect of such infinitesimal distortion of the pencil of rays along the null geodesic from the source 
to the observer. 

Jacobi fields appear naturally in the context ofgravitational lensing They proved particularly useful in 

understanding spacetime singularities (see ||, and ||). More recently, Jacobi fields have been used to understand 
issues in the approach to the Einstein equations referred to as null surface formulation Jl(],|| . The feature of Jacobi 
fields that we focus on in this paper is the fact that the geodesic deviation equations, governing Jacobi fields, acquire a 
particularly simple form in terms of the dynamics (the optical or Sachs equations) for the optical scalars (convergence, 
p, and shear, a) of the null geodesic congruence. 

In this paper we derive a relationship between the distortion of the pencil of rays of the elliptical source and the 
optical scalars of the geodesic congruence. This relationship applies in arbitrary spacetimes, without reference to any 
approximation. Even though our derivation makes use of a non-perturbative definition of image distortion based on 
Jacobi fields, the Sachs equations allow us to eventually eliminate the Jacobi fields and obtain a relationship purely 
between the change in shape of the pencil of rays and the optical scalars. Via the optical equations, this can be 
interpreted as a cause and effect relationship directly between the curvature of the spacetime and the distortion of 
the image. 

Here, as in our companion paper we assume that the source being imaged does not lie across a caustic of the 
past lightcone of the observer. 



II. SUMMARY OF NON-PERTURBATIVE IMAGE DISTORTION 

We summarize, for easy reference, the salient aspects of our approach to image distortion as developed in ||, and 
as represented in Fig. ^ (Readers familiar with Q may want to skip this section). A non-perturbative approach to 
distortion requires, as a fundamental tool, a non-perturbative lens mapping: a mapping from the image location to the 
source location that does not rely on weak-field nor thin-lens regimes. As explained in fll]-|l3[| , such a lens mapping 
can be obtained, in principle, from the expression 

z a = F a (z^r),s,e,cj ) ) (5) 

which gives the coordinates z a of source points on the past lightcone of the observer, located on the worldline z^ (t), 
in terms of the null geodesic that connects the source with the observer. Eq. (|J) can be obtained, for instance, by 
integrating the null geodesic equation, or by solving the eikonal equation Jl4|,|l5| . (Tens equations in generic spacetimes 
without the thin- lens approximation are also being considered by other authors jl6],[l7|]. In particular, see Jl^] for exact 
lensing in Kerr spacetime.) The angles {9, 0) represent the direction of the null geodesic at the observer's location and 
specify the angular location of the image on the celestial sphere in standard spherical coordinates, whereas s gives the 
parameter distance of the source to the observer along the null geodesic (it can be thought of as an aflme parameter) . 
The tangent vectors to the null geodesies in the lightcone are 

dF a 

£ a = -E-- (6) 
os 

Associated with each null geodesic, there is a pair of parallel propagated spacelike vectors, (e",ef), which span the 
space of spacelike vectors orthogonal to £ a , and which allow us to compare angles at two different locations along a 
null ray. The parallel-propagated basis is defined by 

t b V h e\ = l b V b e a 2 = 0, (7) 

e v £ = e 2 -£ = 0, (8) 

e r e x = e 2 - e 2 = -1, (9) 

e r e 2 =0. (10) 
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Even though such parallel propagated basis are only defined up to a fixed rotation, in principle a unique such basis 
could be picked in the case that the electromagnetic radiation emitted by the source is polarized. In such case, the 
polarization vector of the radiation is parallel transported jl) and defines for us one leg, say ef, of our basis. 

The situation that we consider is that of a small source located at the value s* along the geodesic. More specifically, 
s* represents a central point in the intersection of the source's worldtube with the observer's past lightcone. By 
assumption, the source's worldtube does not intersect the caustic of the lightcone, therefore the intersection with 
the lightcone is continuous and produces a single (distorted) image. The source's visible shape is defined by this 
intersection and, as a set, is connected to the observer by a pencil of rays. If the source is small, the pencil of 
rays consists of a bundle of null geodesies that are neighboring to a fixed null ray from the source's "center" to the 
observer. Points on the pencil are thus reached by connecting vectors (Jacobi fields). These are solutions to the 
geodesic deviation equation. Given two linearly independent solutions, (Mf (s), Mf (s)), any other Jacobi field Z a can 
be expressed as 

Z a (s)=aM?(s)+/3M$(s) (11) 

where a and (3 do not depend on s. One natural basis of Jacobi fields is found by taking derivatives of the lens 
mapping, in the form 

~ Bz a ~ 1 Bz a 

oO sin 6 ocp 

Although all Jacobi fields associated with the observer's past lightcone vanish at the observer's location (the apex of 
the lightcone), these two are such that their s— derivatives are orthonormal at the observer's location. As s varies along 

the null geodesic, their components along (e^ef) form an s— dependent matrix that we refer to as J,-(s) = M^e\. 

We see that, by construction, s Jj(s) — * 5* as s — > 0. 

This matrix evaluated at the source's location s* defines for us the Jacobian of the lens mapping, denoted *J. 
Namely, the Jacobian of the lens mapping is 

*J = J(s*), (13) 

and is its inverse *J is our magnification matrix. Notice that, in a weak-field, thin-lens regime, our Jacobian *J is 
related to the Jacobian A of the thin-lens theory via 

*7 

A=-, (14) 

which is due to the fact that we have not scaled our Jacobi fields at the source's location. (In a generic spacetime, 
i.e., in the absence of a flat background, there is no geometric meaning to such a scaling.) 

For ease of describing the source's shape, it is convenient to use a basis of Jacobi fields (Mf,Mf), rather than 
(Mf , Mf ), that is identical to the parallel propagated basis (ef, e%) at the source's location. See Fig. ||. We introduce 
the notation that any quantity defined or evaluated at the source is preceded by a *. We thus have 

*Mf = *e?, (15) 
*Mf = *e a 2 . (16) 

In terms of (Mf , Mf ) we have that 

Mf = (tT^Mf ( s ) (17) 
If we define the components of Mf in the parallel propagated basis ef by 

J){s) = M?(s)4, (18) 
where e J a is the dual to e", we see immediately that 

j}( S *) = (*j _1 )fM;v a = ar 1 )*^; = 5*, (19) 

as required by Eqs. (|l5|)-(|l6|). One also sees, from lim s ^ s_1 ^j( s ) = fy, that 
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hmij}( S ) = (*J (20) 

an important result that allows us to map the source's shape to the image's shape on the observer's celestial sphere. 
This result is fundamental to the construction in this paper, since it allows us to obtain the magnification matrix by 
continuously moving along the null geodesic from the source to the observer. 

It is important to notice that the vectors (e", e 2 ) can be expressed in terms of (Mf , M 2 ), by 

el = a cos(A + v)M% + cos AM 2 a , (21) 
e a 2 = a sin(A + i/)Mf + sin AM 2 a , (22) 

where 



-M 2 -M 2 



Mi ■ MiM 2 • M 2 - (Mi • M 2 ) 2 



(23) 
(24) 

Mi • MlM 2 • M 2 - (Mi • M 2 ) 2 ^ ; 

Mi-M 2 

cos t/ = ; (25) 
VMi-MiM 2 -M 2 



-Mi -Mi 



and A is determined, up to a constant, by 



U , Ml ' Ma j-logf^l (26) 



2i/Mi • MiM 2 • M 2 - (Mi • M 2 ) 2 ^ V M i' M 2 
At the source, by Eqs. (p^|)-(^6|), we have 

*a = l, *0=1, * v =~\ * X =\ ( 2? ) 

A small elliptical source at s* of semiaxes L + and i_, oriented so that the semimajor axis lies at an angle 5* with 
*el, is described parametrically by 

*Y a {t) = Y\t) *e c l + Y 2 {t) *e a 2 = y x (i) *M? + Y 2 {t) *M 2 a , (28) 

where the components Y l (t) (for t € [0, 27r]) can be specified as 

y 1 = cos i cos*<5 - sin t sin*<5) , (29) 

Y 2 = (*£ cos t sin*<5 + sin t cos*<5) . (30) 

Here 1Z is the ratio of the semiaxes (1Z = L+/L-) and A is referred to as the area of the ellipse, even though it is 
defined simply as the product of the semiaxes (^4 = L+L_). We can follow the lightrays that connect each point of 
this source with the observer by defining the connecting vector Z a (s,t) as 

Z a (s,t) = Y 1 {t)M1{s) + Y 2 (t)M%(s) (31) 

The image is thus obtained by projecting Z a (s, t) along the directions (e", e 2 ), in the limit as s — > 0: 

X*(t) = lim ly(f)M»4 = (*J _1 )}^(t) (32) 

In order to study the distortion of the image, compared to the shape of the source, we define the shape parameters 
to be the area of the ellipse A(s), the ratio of its semiaxes TZ(s) and its orientation S(s), for each s as the ellipse is 
carried by the null geodesies towards the observer. Defining, for short, the following 
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a(s) = *TZ 2 ((cos*5) 2 Mi • Mi + (sin*<5) 2 M 2 • M 2 + 2 sin*6cos *SM 1 ■ M 2 ) , (33) 
b{s) = -(1 - *ft 2 )(M r Mi + M 2 -M 2 ) + 1(1 + *K 2 )cos(2*<5)(M 2 - M 2 - M r Mi) 

-(l + *ft 2 )sin(2*<5)Mi-M 2 , (34) 

c(s) = *ft(sin(2T)(M 2 - M 2 - Mi- Mi) + 2 cos(2*£)Mi • M 2 V (35) 

it can be seen Q that the shape parameters can be obtained from the Jacobi vectors (Mf , M 2 ) as 

.4 = ^(Mi • MiM 2 • M 2 - (Mi • M 2 ) 2 ) V2 (36) 

K={ M --fHX'" (37) 

acos(A + v){*TZcofit + sin*5 + sini + cos*8) — /3cos A(*7?.cos£ + cos*<5 — sini + sin*5) 
/3sin A(*7?. cosi + cos*i5 — sin^ sin*(5) — a sin(A + ^)(*7?.cost + sin*<5 + sini + cos*<5) 



where t + is the value of t that corresponds to the major semiaxis of the ellipse, which is given by 

1 \ 




(39) 



[Other authors [[l9| have studied the rotation of the image of an elliptical source, corresponding to our Eq. j3q ) , but 
in a different physical context (the appearance of the polarization of radio emission by sources at high redshift).] 
It is useful to define, as well, the angular analog of the area of the ellipse, namely, the solid angle of the pencil of rays: 

Q(s) = 4 (40) 

The image subtends a solid angle f2/, and its semiaxes have a ratio IZj and orientation Si given by 

A(s) 

fi/ = lim — Y- = lim ills), (41) 

s— >0 s s^O 

Hi = K{Q), (42) 
Si = 6(0). (43) 

The distortion of the image with respect to the source is measured by the distortion parameters 25q, and T)$ 
defined as 

s*\ 4 A 2 tl 2 



0nS &UJ w-^Hw- 1 * (44) 

D a = 5(0) - *S. (46) 

If any of the distortion parameters is non vanishing, the image is said to be distorted with respect to the source. The 
reader is referred to ||] for complete details of the construction in this Section. 

As defined, the distortion parameters are expressed in terms of the Jacobi basis (Mf , Mf). Even though it might be 
tempting to colloquially refer to Dfj as the convergence, and to Dr as the shear of the image, here we are interested in 
obtaining a relationship between the distortion parameters and the proper convergence p and shear a of the lightcone. 
In the remainder of this paper we show that the dependence on the Jacobi fields can, in fact, be eliminated in favor 
of the convergence p and shear a of the lightcone. 

Before we demonstrate the relationship between the distortion parameters and the optical scalars in full generality, 
we take a preliminary step. We develop the relationship in the case of a thin lens treated non-perturbatively. This 
amounts to applying the formalism as developed in this section to the case that the curvature of the spacetime along 
one light ray is a Dirac-delta function. 
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III. THIN LENSES TREATED NON-PERTURBATIVELY 



In this Section, we assume that the spacetime is given, and that it is flat everywhere in the region of interest except 
at a surface, that will be interpreted as a lens "plane" . We will not be concerned with the geometry of this surface, 
which in fact may or may not be a plane. We actually restrict attention to the spacetime in the neighborhood of a 
given lightray that connects a small source to the observer, and assume that the situation is generic to the observer's 
lightcone, except in the neighborhood of a caustic. This situation can be modeled with a distributional curvature 
tensor along the given lightray, with support at the lens plane. 

Jacobi fields Z a are solutions to the geodesic deviation equation 



£ c V c £ b S7 b Z a = R\ cd l b Z c t d 



(47) 



along a fixed null ray with tangent £ a . This is an ordinary second-order differential equation for the components of 
the Jacobi field, where the curvature of the spacetime along the ray is considered as given. Because the space of 
solutions is two-dimensional (in the case of null geodesies), then this is essentially a problem of two equations for two 
components. In the following, we summarize a standard choice |||7]] for the reduction of (47) to a linear problem in 
two dimensions. 

Firstly, a parallel propagated basis transverse to the null geodesic congruence is given. This is the basis (e", ef) 
introduced in the previous Section. However, it is customary to "complexify" the basis, mainly for easier handling 
and notation. We define 



b — , — 

V2 



(48) 



which consequently satisfies m-m = —1 and ra-l — m-m = £ b Vbm a = 0. We consider (m a ,m a ) as our parallel 
propagated basis. Similarly, we arrange our basis of Jacobi fields (Mf, M% ) into a complex basis 



M a 



7T~ 



which can now be expressed in terms of (m a ,m a ) by means of two complex components (£,??) via 

M a = £m a + rjm a 



(49) 



(50) 



(Here, £ and rj should not be confused with the source location on the source plane and the image's location on the 
lens plane, which, as explained in Section I, is the notation used in JjJ], unfortunately). Substituting Eq. ( |5C| ) into 
Eq. (4?]) and contracting with rh a we obtain 



^ 2 C = -(« + *o?y), 

where D = £ a V a , $oo = \R a b£ a £ b and = C a bcd a m b £ c m d . Similarly, contracting with m a yields 

We can arrange these two equations into matrix form as 

D 2 X = QX, 

where X is a matrix containing the components of the connecting vectors 



X 



and Q is a matrix containing the relevant components of the curvature: 



Q = 



$00 *0 
*0 $00 



(51) 
(52) 
(53) 

(54) 
(55) 



Equation (|5^), which is equivalent to Eq. (^7) 
from given free initial values, rather than by 



represents an alternative way to obtain Jacobi fields, by integration 
derivation of the lightcone. However, there is still an alternative version 



G 



of (|53|) for Jacobi fields which is a first-order formulation involving the optical scalars, rather than the curvature. To 
obtain this formulation, we notice that all connecting vectors satisfy 

i a V a Z b = Z a W a £ b , (56) 



namely, they are Lie-dragged along the congruence. Putting Eq. fl50| ) into Eq. (|56|) and contracting alternatively with 
m a and rh a we obtain 

DZ = -(j>S + aTi) (57) 
D V = -(a£ + m) (58) 

with 

p = m a m b V b l a (59) 
a = m a m b V b l a (60) 

which is equivalent to the matrix equation 

DX = PX (61) 

with 



P ° 
a p 



(62) 



Eq. (J6l]) is a first-order differential equation for the Jacobi fields in terms of the optical scalars: the convergence p, 
and the shear er, of the geodesic congruence t a . As for P itself, it can be obtained directly from the curvature by 
virtue of the Sachs equations for the optical scalars, namely, 

DP = Q + P 2 , (63) 

which is obtained by taking a derivative of Eq. ( |6l|) and using Eqs. (H) and © to eliminate X from the resulting 
equation. Thus we think of the optical scalars p and a as given, for the purposes of solving for the Jacobi fields in 
Eq. @. 

We now specialize our discussion to the particular case of a thin lens. Namely, the curvature Q is zero outside of 
a single, two-dimensional, spatial surface. Under this assumption, there is only one value of s along all relevant null 
rays in the past lightcone of the observer where the value of the curvature is non-zero, i.e., the point where the null 
geodesic intersects the lens "plane" . This value will depend on the choice of (6, </>) labeling the null geodesic but we 
will simply call this value s — L < 0. For a given ray, the curvature matrix may be written as 

where ($oo, |*o|i <p) are rea l constants. 

In the spirit of our current program, the basis (Mf , Mf ) that we wish to solve for lies on the past lightcone of the 
observer and coincides with (e", ef ) at s*. Translating these two conditions into conditions for X, by Eqs. (J5fj), ( |49| ) 
and (f4§|), we obtain two boundary conditions on the solutions to Eq. (|53|), respectively of the form 

X(0) = (65) 
X(s*) = I (66) 

where I is the identity matrix. The vanishing boundary condition at the observer ensures that the observer receives 
all the neighboring lightrays. The unit boundary condition at the source ensures that the two Jacobi fields that we are 
solving for are orthonormal and aligned with the basis (e°, e§) at the source location. [It is worthwhile noticing that 
we could have, just as well, imposed boundary conditions that would allow us to obtain the basis (Mf,Mf ). These 
would require that X be vanishing and DX be the identity at s = 0. By Eq. ( |l3| ) and Eq. (|2(]), both procedures are 
equivalent for the purposes of finding the magnification matrix. ] 

At all points but s = L, Eq. ([33]) reduces to D 2 X = 0, thus the solution has the form of two linear functions of 
s at both sides of the lens plane. We denote these functions as Xo and X$ for the observer's side and the source's 
side, respectively: 
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J X = sK L<s<0 
A ~ \X s = (s-s*)N + 1 s*<s<L ^'1 

where K and TV are two constant matrices which are determined by matching conditions on the lens plane. The 
matching conditions, according to Eqs. (|5^) and (pi}), are continuity of X and a jump in the first derivative of X, 
namely 

X (L)=X S (L), (68) 
DXo(L) =DX S (L)-TX S (L). (69) 

This allows us to obtain K from the curvature matrix T as 

K = (s*I - L(L - s*)T] 1 (70) 

and N from K and T as 

N = K + LTK (71) 

Our goal is, however, to find K in terms of the optical scalars, so we need to express T in terms of P. This can be 
done very quickly here because of the thin-lens regime that we are imposing. From Eq. ( |6l|) we obtain K = —sPoK, 
or 

Po = -- (72) 
s 

on the observer's side or the lens plane, as it should be since on the observer's side the lightcone is the same as in flat 
spacetime. On the source's side, however, we obtain 

N = -P s ( K (s-s*)N + l) > (73) 



which is valid at all values of s on the source's side. In particular, we evaluate it at s — L and use Eq. ( p8|) on the 
right-hand-side and Eq. (Fry) on the left-hand-side, obtaining 



K + LTK = -LP S (L)K, (74) 

or 

Ps(L) = ~-T. (75) 
This means that the curvature is related to the values of the optical scalars at the lens plane (on the source's side) by 

/>s(£)=- (jr + 5oo) , crs(L) = -|*o|e lv . (76) 

Substituting Eq. @ into Eq. @ we obtain 

k „- , )PsW y> tfr-Z. r + ( ^*& ) ™ 

with 

A= (l + (L-s*)p L ) 2 ~(L~s*) 2 <T L a L , PL = ps{L) : a L = a s (L). (78) 

For completeness, we display the form of the shear and divergence in the source's side of the lens, which can be 
obtained directly from Eq. (p3) by setting Q = (or, less directly, from Eqs. (|7l|), (|73|) and (|77|) we can solve for Ps)- 



(pi - °lol)\Pl - (s- L)(p 2 L - a L a L )) 

Ps(s) = — ^ - 2 >- (79) 

{PL - (s - L)(p 2 L - a L a L )J 



os{s) = (Pl-^)a L (8Q) 



(p L - (s- L){p\- <J L a L )^j 
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We can now read off the complex components of the complex connecting vectors M a and use them to obtain the 
real connecting vectors M-f and M%. We have 

M a = (j + i+l+j) e q + (R^A_A±H e a (81) 

Mi = (AdL+l^l e a + tf±izlZ^^ (82) 
2i 2 

which, with Eq. @ and Eq, @, yields 

Ml = ^(1 + (I - »*)Pi — (£ — cospK - JL((L - s*)\a L \ sm^K (83) 

M 2 ° = -^((i - **)KI sin^e? + -^.(1 + (I - S *) P£ + (I — Okil cos^)e§. (84) 



Consequently, by Eq. (18), the components of the matrix J on the observer's side are 



s 



1 + (I - s*)p L — (I — s*)|cr L | cosp -(I - s*)ki| sinp 



J= LA I (85) 

LA 1 -(L-s*)|<7 L |sin<p l+(I-s*)p L + (I-s*)|a L |cos^ 

By Eq. (p0[), our magnification matrix *J , which maps the source's shape intoO the image's shape, is thus 

, i / l + (L-s*)p L - (L- s*)\a L \ cos (p -{L - s*)\a L \ simp 

J 



LA 



-(L - s*)\a L \ sinp 1 + (I - s*),o L + (I - s*)|ctl| cosp 



To obtain the Jacobian matrix A of the thin-lens approach, we use Eq. ((Tj), which takes care of the scaling of the 
source vector. Inverting Eq. ( |S6| ) and dividing by s* we obtain 

T ( 1 + (L-s*)p L + (L-s*)\a L \ cos <p (I - s*)\a L \ simp 

A =-A I- (87) 

(I - s*)|cr L |sinp 1 + (L-s*)p L - (I - s*)|cr L |cos<p 

By comparison with Eq. (^|) , we can make the following identifications 

«= ~* L \ l + Lp L ), (88) 

71 = ^(s* -L)\a L \ cos cp, (89) 

72 = 4(s*-I)k i |sin^. (90) 



In the weak-field, thin-lens approach to lensing, k and 7 = Vo+tI are referred to as the convergence and shear. 
We can now see how they relate to the values of the optical scalars of the lightcone - the convergence and shear-, 
evaluated at the lens plane on the source's side. We can verify, by inspection, that k, ji and 72 vanish in the case 
of flat space, where pl = — 1/1 and ul = 0, which means that there is no distortion in the absence of a lens (in 
the thin-lens regime, it makes sense to define distortion away from the distortion associated with flat space, which is, 
naturally, the convergence of the lightcone) . 

The generic case in which the curvature along a lightray is continuous also offers a relationship between image 
distortion and the optical scalars p and a, certainly not in terms of the values of the optical scalars at particular 
points. We phrase this relationship in the context of the shape parameters, rather than the magnification matrix. 



IV. GENERIC NON-PERTURBATIVE APPROACH TO IMAGE DISTORTION 

In this section we generalize the results of the previous section to the case of generic spacetimes. Our goal is to find 
a direct relationship between the optical scalars and the distortion of the cross-section of the pencil of rays connecting 







the observer to an elliptical source. Such a relationship will be interpreted as the direct effect of spacetime curvature 
in the distortion of small elliptical images. 

That such a relationship exists can be seen from the following. By Eqs. (^5|)-(p6|), the distortion parameters Dq, T>tz 
and Dg are expressed in terms of the scalar products of the Jacobi basis (M", M%). The Jacobi basis vectors (Mf , M% ) 
along a fixed lightray need not be derived from the lightcone, but can be obtained directly from the optical scalars 
by means of Eq. (|6l|), as we showed in the previous section for a special case. Therefore, the distortion parameters 
can be thought of as functionals of the optical scalars. Here we present one way of obtaining an explicit relationship 
between thedistortion parameters and the optical scalars. The relationship that we obtain is a system of first-order 
ordinary differential equations for the evolution of the shape parameters along a fixed lightray, where the optical 
scalars enter as known sources. This system of equations can be thought of as an extension of the focusing equation 
for the determinant of the Jacobian matrix in lensing [jlj. Significant differences are, however, that the equations are 
of first order (instead of second order) and that all three parameters are evolved, not just the area of the beam. 

To obtain this result we start by taking a derivative of the shape parameters with respect to the affine parameter s 
and evaluating the derivative at the value s*, at which the Jacobi basis is orthonormal. This can be regarded as an 
infinitesimal distortion, in a sense to be made specific. We have, from Eq. (^6|), for the infinitesimal distortion of the 
area: 



dA 2 



ds 



-A 2 -(M r M 1 +M 2 - M 2 ) 
ds 



The infinitesimal distortion of the semiaxes ratio can be obtained from Eq. p?n by 



dTZ 2 



ds 



dTZ 2 da dTZ 2 db dTZ 2 dc 



da ds db ds 



dc ds 



where (a, 6, c) are the functions of s defined by Eqs. (p3|)-(p5|). From Eqs. 

a(s*) = -*K 2 , 
b(s*) = *TZ 2 -l, 
c(s*) = 0. 



|), ( p4[) and (|35| ) we can see that 



By calculating the partial derivatives of Eq. (|37j) and using Eqs. (| 

on 2 



|) we obtain 



da 
dTZ 2 



db 
dTZ 2 



1Z Z 



*TZ 2 , 



1, 



dc 



(91) 



(92) 



(93) 
(94) 
(95) 



(96) 
(97) 
(98) 



Taking an s— derivative of Eqs. (|3 
@, in Eq. H) yields 



5[), evaluating it at s* and using the resulting expression, together with Eqs. (t 



dTZ 2 



ds 



= 1Z 2 ( cos(2<5)-^(A/ 2 - M 2 - Ah- Ah) - 2sm(25)-^M 1 - Al 2 

» \ ds ds 



(99) 



The infinitesimal distortion of the orientation parameter is more complicated but can be obtained in a straightforward 
manner via a similar procedure: 



d5_ 
ds 



ds 



(100) 



Notice that Eqs. (^lj), ( p9| ) and ( |100[ ) depend on the derivatives of the scalar products between our Jacobi basis 
vectors. In order to eliminate these in terms of the optical scalars we find first an expression for the optical scalars 
in terms of them. These can be obtained directly from their definitions, Eqs. ( |5S| ) and (|60|), using Eq. ( fi"8"| ) and 
Eqs. ©-(IH). We have 



m a = ^- e l{XJrv) All + 4?e lA M 2 a 
V2 



V2 



(101) 



10 



1 

2 

a=- [ o 



p=- [a^M^Ml + fi z M%M b 2 + 2a(3 cos v{M%M% + M{M%)) V a h 

2 e 2i(\+v) M a M b + 0i e 2iX M a M b + a p e W>>+») (M? M$ + MlM%))V a tb 



(102) 
(103) 



where we have used the fact that Vf,£ = V a ^t- In order to put Eqs. ( PL02I) - ( pL03|) in a convenient form, we use the fact 
that that the connecting vectors M-f , M 2 are Lie-dragged by the geodesic vector £ a , namely, Eq. (|5^) applied to both 
Mf and M 2 a : 



(104) 
(105) 



This allows us to trade the covariant derivatives of £ a for ordinary derivatives of the scalar products between the two 
connecting vectors: 



1 



P=-A a* — (M 1 -M 1 )+f3 z — (M 2 -M 2 )+2a/3cosv—(M 1 -M 2 ) 



ds 



ds 



ds 



J2i\ 



a = — cte Mv — (M r M 1 ) + /3 2 — {M 2 -M 2 ) + 2a(3e lv — (M r M 2 ) 



ds 



(106) 
(107) 



For c ompl etene ss, w e display the final expressions for the optical scalars in terms of the connecting vectors, alt houg h 
Eqs. (106) and (107) are sufficient for the purposes of this section. Substituting Eqs. (|23|), ( pi| ) and (25) in Eq. (106), 
the expression for p becomes 



p = - 7 — log(Mi • Mi M 2 ■ M 2 - (M r M 2 f) 
4 ds 



(108) 



Using Eqs. (|23[), ( |24| ) and ( J25| ) in Eq. (107), we obtain an expression for the magnitude of g: 

(M 1 -M 1 M 2 -M 2 f 



GG 



4(Mi • M X M 2 ■ M 2 - (Mi • M 2 ) 2 ) 2 



log- 



Mi -Mi 
M 2 -M 2 



4(Mi-M 2 ) 2 



log 



Mi • MiM 2 • M 2 ds V M • Mi 7 ds 



Mi-M 2 



log 



M r M 2 
M 2 -M 2 



(109) 



In order to obtain an expression for the phase <p of the shear (ct = |cr|e iV ), we first calculate the phase of g = |cr|e l< ' 5 , 
substituting Eqs. (||), (||) and (||) into Eq. (|l07| ): 



tan<£ = 



2M r M 2v ^f72^1og(f^) 



MiMiM 2 M 2 f l 0g (Mf±h\ + 2(M 1 -M 2 )^ log (f^) 



The phase of the shear is thus 



= 2A + ^, 



(110) 



(111) 



where A is determined by Eq. (|26|). The set of equations Eqs. (108), ( |l09| ) and (111) is the equivalent of Eq. ( pl|) o f 
the previ ous Section. For our pur pose s, however , we use the preliminary form of the equations given by Eqs. (106) 
and (107). We evaluate both Eq. (106) and Eq. (107) at s* , using Eq. (|27| ) for a,/3, v and A. We obtain 



P* = 7 -j- (Mi -Mi +M 2 -M 2 ) 
4 ds y ' 

Id, , 

a-* = (Mi -Mi - M 2 -M 2 + 2jM r M 2 ) 

4 ds v " ' 



(112) 
(113) 



Subtituting Eq. (|TT|) and Eq. ( pT3| ) into Eq. ©,Eq. © and Eq. ( |l00| ), we obtain 
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dA 2 
ds 
dTZ 2 



ds 
d5_ 
ds 



= -4*p % 2 

= -2(cos(2*<5)(W - isin(2*<5)(*cr- *a))*1Z 2 
{1Z 2 + 1) 



2(*K 2 - 1) 



^ ( sin(2*<5) + V) + % cos(2*<5) (*a - &)) . 



(114) 
(115) 
(116) 



Since these equations hold for arbitrary fixed value of s*, we can drop the stars and consider them as differential 
equations for the shape parameters 



dA 2 
ds 

dTZ 2 
ds 



= -Ap A 2 

= -2(cos(2<T)(cr + a) - ism{25){a - a))n 2 



dS _ (K 2 + l 
ds ~ 2(1Z 



j—t ( sin(2<5)(<7 + a) + i cos(2<5)(ct - ajj . 



Notice that, as a consequence of Eq. (117), the solid angle SI = A/s 2 satisfies the differential equation 



dn2 J r \,v 



(117) 
(118) 
(119) 

(120) 



which shows that in flat space, where p — — 1/s, there is no change in the solid angle of the image 0/ as compared 
to the solid angle of th e sourc e *f2 , since there is no change in the solid angle of the pencil of rays Jl(s). We consider 
alternatively Eq. (117) or Eq. (120) as equally relevant in order to refer to infinitesimal area distortion or infinitesimal 
solid angle distortion. However, it is the infinitesimal solid angle distortion that is relevant to total distortion, as we 
see immediately below. 



Equations (117) (or equivalently (120)), (118) and (119) constitute our desired direct relationship between image 



distortion and the optical scalars. In these equations, the optical scalars p and a are thought of as given, since they 
are obtained, along a fixed null ray, from Eq. (|6^) where the spacetime curvature acts as the source. As expected, 
the convergence p (also referred to as divergence and expansion elsewhere) enters only in the equation for the area 
distortion. Notice, however, the role of both the real a nd im aginary parts of the shear a in both the orientation 
and semiaxes ratio distortions. Notice, as well, that Eq. (119) breaks down when the ratio of the semiaxes takes the 
value 1, which reflects the fact that the orientation is ill defined at the points where the shape of the pencil of rays is 
circular. 

The finite distortion parameters ^iz and are obtained by direct integration of the infinitesimal distortions: 



1 f° dn 2 , /f2 z 



i 



1 f° dK 2 , 
Stc = I ^^ds 



*R? J s » ds 

Ds = [ ^-ds = (5/ - *<5. 
ds 



(121) 
(122) 
(123) 



As discussed in [0], we say that there has been image distortion if any of the three dist ortion parameters is non-zero. 
As we noted, in flat space, p — — 1/s, so there is no change in the solid angle in Eq. (120) an d the distor tion parameter 
£>si is zero. Likewise, a is zero in flat space, so D-r, an d D s are also zero by virtue of Eqs. (118) and (119). 

In order to illustrate the significance of Eqs. (117), (118) and ( |119| ), in the following subsection we specialize to the 
case of a Schwarzschild spacetime. 



The Schwarzschild case 



Consider a Schwarzschild spacetime in standard Schwarzschild coordinates 



ds 2 



2m i 
1 ) dt 2 

r 



2m 



dr 2 - r 2 (d6 2 +sm z 0d(j) 2 ). 



(124) 
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The deflector and the observer define a coordinate line that is fixed for the lensing problem, referred to as the optical 
axis. For simplicity, we choose this to be the z— axis. Because the spacetime is spherically symmetric, the lensing 
problem (defined essentially by the past lightcone of the observer located on the z— axis) has axial symmetry around 
the optical axis. Every geodesic emitted "backwards" in time from the observer remains on a coordinate plane with 
a fixed value of <f>. Therefore, there is a connecting vector, Mf , that lies on the plane <f> = constant (i.e, has no 
component along d/d(f>) and there is a second connecting vector, M$, that joins geodesie s in adjacent <p = constant 
planes (i.e., can be taken to have only a component along d/d<j)). With the metric, Eq. (124), these two connecting 
vectors are orthogonal everywhere on the lightcone. In terms of these connecting vectors we have that 



1 d / Mi' Mi 
Ids ° g \M 2 -M 2 



which follows from Eq. (107) by using the facts that v = constant 

and (§( 



applied to Eqs. 



(125) 

= 7r/2 and A = constant — (by choice), both of 
h. Thus the shear is real (up to a constant phase 



which are consequences of Mi ■ M 2 
that has been chosen as zero). 

If the shear is real, the differential equation for the orientation parameter admits the solution S(s) = 0. Thus, if 5 
starts out with the value zero, it remains zero along the given null ray. With S — 0, the equation for the ratio of the 
semiaxes, Eq. (|118|), reads 



dK 2 
ds 



(126) 



with solution 



K 2 = *TZ 2 exp -4 



(127) 



If one uses Eq. (125), it can be seen that this agrees with the direct calculation of the ratio carried out in our 
companion paper m. 



On the other hand, if 5 is not vanishing at the start, then it does not remain constant, by Eq. (HE), and we have 
a system of coupled differential equations for 1Z and 8. This reflects the fact that the images of sources that respect 
the axial symmetry of the lensing problem also respect such symmetry However, the images of sources that do not 
have axial symmetry around the optical axis will be rotated, as described in our companion paper . 



V. CONCLUDING REMARKS 



We have attempted to provide a unified description of image distortion in a non-perturbative (exact) manner within 
the framework of general relativity. We have tried to formalize our description in a manner as close as possible to the 
standard (weak-field thin-lens) treatment with the hope that it will be useful in deepening the current understanding 
of the gravitational lensing phenomenon. 

We have related the convergence and shear used in standard lensing to the actual values of the geometrically defined 
optical scalars (convergence and shear) at the lens plane. 

The roles of our non-perturbative distortion parameters can be compared with those of the convergence and shear 
from the thin-lens approximation. Both descriptions of image distortion are defined such that they vanish in flat 
space, so that a non-zero value indicates a deviation from an "unlensed" image. The k in the thin-lens app roximati on 
is analogous to 2)q in that they both give an indication of an overall change in angular size. The 7 = vTi+tI m 
Eqs. ( p0[ ) indicates stretching and squeezing of the image along perpendicular axes, as D-ji does. Lastly, the ratio 
72/71 gives the orientation of the axes of the image, as 1>s does. The analogy is not perfect, but the issue is that 
there are as many parameters to completely describe image distortion in one case as in the other. 

Clearly, image distortion is ultimately governed by the geodesic deviation equations for the connecting vectors of 
the pencil of rays between the source and the observer. This is a system of second-order equations for four variables 
(the two components of two linearly independent connecting vector fields which span the space of solutions, or, 
equivalently, the four components of the 2-dimensional linear map D\ that maps the initial values of a connecting 
vector Zj(0) into its values at a point s via Zi(s) = D J i Zj(Q)). Still, the distortion of elliptical images has only three 
relevant parame ters . Defin ing t he th ree distortion parameters and using their evolution equations in terms of the 
optical scalars, (118), (|llSj ) and ( |l20|) , allows us two advantages: 1) the problem is a system of equations for the right 
number of relevant quantities, and 2) the system is of first order, rather than second. 
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In the end, defining the distortion parameters allows us to reduce, in a way, the geodesic deviation equations down to 
the specific problem of elliptical images. We feel that they may become useful in studies of weak lensing with multiple 
lens planes, cosmic shear [^0| and distortion of high-redshift sources [^lj. In particular, meaningful observables can 
actually be defined out of the shape and distortion parameters. Work on this issue is in progress and will be reported 
elsewhere. 
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o 




FIG. 1. The lensing problem in standard gravitational lensing. The point O represents the observer. The point S represents 
the source. The lens mapping can be interpreted as a map that takes the point / in the lens plane into the point S in the 
source plane. 
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observer' s\ 
worldline 



source' s 
worldtube 




FIG. 2. Diagram of image distortion. The source's worltube intersects the observer's past lightcone in a region free of 
caustics. The source's visible shape is denned by the intersection. The pencil of rays between the source's shape and the 
observer carries the shape of the source into the shape of the image, on the observer's celestial sphere. The pencil of rays can 
be described by geodesic deviation vectors (connecting vectors of the observer's lightcone) from a central null ray connecting 
the center of the source's shape to the observer. 
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FIG. 3. Vectors associated with the null geodesic that joins the source to the observer. The vector i is tangent to the null 
geodesic. The vectors (ei,e2) are spacelike, orthogonal to £, orthonormal and parallel progagated along the geodesies. The 
vectors (M\,M2) are linearly independent Jacobi fields which coincide with (ei,e2) at the location of the source (not shown). 
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